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Abstract. Circuit models involving bulk currents and inductances are often used 
to estimate the energies of coronal magnetic field configurations, in particular con- 
figurations associated with solar flares. The accuracy of circuit models is tested 
by comparing calculated energies of linear force-free fields with specified boundary 
conditions with corresponding circuit estimates. The circuit models are found to pro- 
vide reasonable (order of magnitude) estimates for the energies of the non-potential 
components of the fields, and to reproduce observed functional dependences of the 
energies. However, substantial departure from the circuit estimates is observed for 
large values of the force-free parameter, and this is attributed to the influence of the 
non-potential component of the field on the path taken by the current. 



1. Introduction 

Solar flares are understood to involve the conversion of magnetic en- 
ergy in the solar corona into the energy of particles and radiation. In 
particular flare energy is believed to derive from non-potential coronal 
magnetic fields, i.e. that part of the magnetic field produced by currents 
flowing in the solar corona (e.g. Tandberg-Hanssen and Emslie, 1988). 

The origin of flare energy in coronal currents has suggested to many 
authors the use of a circuit analogy (e.g. Alfven and Carlqvist 1967; 
Spicer, 1982; Somov 1992; Melrose, 1995; Melrose 1997; Longcope and 
Noonan 2000; Khodachenko and Zaitsev, 1998; Khodachenko, Haeren- 
del and Rucker, 2003). It is well known that the energy of a system 
of N distinct current-carrying circuits may be written (e.g. Jackson, 
1998) 

N N 

£ = oE^ + EE Mijijj, (i) 



where 



=i j>i 



H ~ 47T/2 J Cl % } Ci 1 |X. 

is the self-inductance and 
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is the mutual inductance, and where J(x) is the current density at a 
point x in space, and C{ is the volume of the i th circuit. Although 
Equations (^Q) — © are quite general, their utility relies on having a 
small number of well-defined circuits with simple geometries for which 
analytic expressions for Li and M^ may be evaluated, and in particular 
for which Li and My are independent of ij and Ij. This is the situation 
for currents flowing in a small number of "wire" loops, i.e. loops with 
a fixed geometry. Circuit models for coronal magnetic fields adopt wire 
loop configurations. For example, in the flare model of Melrose (1997) 
based on reconnection between current-carrying magnetic flux tubes, 
each flux tube is represented by fixed toroidal structures carrying a 
uniform axial current. The self-inductance of such a structure is (e.g. 
Landau and Lifshitz, 1960) 

Li = noaiCi, Q = In — - - -, (4) 

n 4 

where Oj is the major radius of the toroid and ri is the minor radius. 
Melrose (1997) also adopts the following approximate expression for 
the mutual inductance between two loops whose centres are separated 
by a distance dij and whose axes subtend an angle dif 

= fi ^^ (5) 

(ai + a,-) 2 + 4 

This formula is not exact but is based on interpolation between known 
results. An exact formula which holds for two axially-aligned line- 
current loops separated by a distance dij is (Stratton, 1941) 



Mij = no(aidj 

where 



d/ 2 




K(hj) - ^E{hj) 



kij 



(6) 



2 Aajaj 

K V {a . + aj) 2 + d 2.> V) 

and where K(x) and E(x) are complete elliptic integrals. 

Real coronal currents do not flow along simple wire loops — they 
are distributed in the coronal volume, and describe complicated paths. 
A common model for current-carrying coronal fields is the force-free 
model, in which the magnetic field B satisfies J x B = 0, (as well as 
V • B = 0) where J = x B. In other words the current density is 
everywhere parallel to the magnetic field. Unfortunately the force-free 
field equation is in general nonlinear and hence difficult to solve. A 
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restricted set of force-free fields — linear force free fields — are easy 
to calculate (e.g. Nakagawa and Raadu 1972; Chiu and Hilton 1977; 
Alissandrakis, 1981; Gary, 1989). For these fields V x B = aB where 
a is a constant, the force- free parameter. Linear force- free models for 
coronal magnetic fields permit detailed descriptions of current paths 
in the corona, although they have limitations as large scale models for 
solar fields (e.g. Sturrock, 1994). 

In this paper the accuracies of simple circuit models for coronal 
magnetic field energies are tested by application to linear force-free 
field configurations. To our knowledge this is the first quantitative test 
of the accuracy of the circuit model in the solar context. For simplicity, 
the comparison is limited to field configurations representing one and 
two magnetic loops. The order of presentation is as follows. In § 2 the 
method is outlined, including the adopted boundary conditions (§ 2.1), 
and the details of the linear force-free (§ 2.2) and circuit (§ 2.3) models. 
In § 3 the results of the calculations for one loop (§ 3.1) and two loops 
(§ 3.2) are presented. Finally in § 4 the results are discussed. 



2. Method 

2.1. Boundary conditions 

The plane z = is taken to represent the photosphere, and the half 
space above this plane (z > 0) represents the coronal volume. We 
consider simple boundary field configurations suggestive of one coronal 
loop: 



B z (x,y,0) = Biexp 
and two coronal loops: 
B z (x,y,0) = Biexp 

+ i? 2 exp 



(x~x 1+ )- 
2a? 



(X - Xih 



2o\ 

(x ~ x 2 ^ 
2a\ 



B\ exp 



— B\ exp 

- B 2 exp 



(x-xi_; 

2a? 



(8) 



(x - Xi- 



2a\ 



(x - X 2 - 



2a\ 



(9) 



respectively, where x = xx + yy. The positions xi± and x 2 ± define the 
loop footpoints. 

It is convenient to adopt non-dimensional units in which lengths are 
expressed in terms of a basic length L, and magnetic fields are expressed 
in terms of the peak vertical field B\ at the base of loop one, e.g. 



Lx, B z = B\B Z , etc., 



(10) 
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where bars indicate dimensionless versions of quantities. 
2.2. Linear force-free model 

A solution to the linear force-free equations is provided by (e.g. Alis- 
sandrakis, 1981) 

r i \ —i(uk — va) i-v.o/ \ 
Uu,v,z) = 2n \ u2 + v2) e kz b° z (u,v), 

i , -\ -i(vk + ua) _fc-=. 0/ , 
by( u ^,z) = Mu2+v2) e kz b° z (u,v), 

b z (u,v,z) = e- k *b° z (u,v), (11) 

where bi(u,v,z) denotes the Fourier transform of B{(x,y,z) in x and y 
(u and v are the wavenumbers corresponding to x and y respectively), 
b®(u,v) = b z (u,v,0) is the Fourier transform of B z (x,y,0), a = aL is 
the force-free parameter, and 



k 



4 vr 2( n 2 +u 2 ) _ s 2j 1 /2_ 



Equation (fTT|) requires (u 2 + v 2 ) 1 / 2 > |a|/(27r). Since large u and v 
correspond to variation on small spatial scales, this solution is termed 
the "small scale solution" to the linear force-free field equations. Of 
course, inversion of this solution in Fourier space involves integration 
over all u and v. The small scale solution only provides a solution 
matching given boundary conditions B z (x,y, 0) if b° z (u, v) is identically 
zero for (u 2 + v 2 ) 1 / 2 < |a|/(27r). In general this is not true, and in 
particular it will not be true for the boundary conditions (JHJ) and 
for any given [although b° z (u, v) is expected to be small for small 

u and v, and identically zero for u = v = provided the net flux is zero]. 
In the following we effectively calculate b^(u, v) for the given boundary 
conditions and then impose b z (u,v) = for (u 2 + v 2 ) 1 ^ 2 < |a|/(2-7r). 

The physical justification for considering only the small scale so- 
lution is that it provides a linear force-free field with finite energy. 
Following Alissandrakis (1981) and Gary (1989), using Parseval's theo- 
rem the energy of the linear force-free field 1)11(1 in the half space z > 
is 

,oo roc \b° z (u,v)\ 2 



Ef = I du i dv u 2 ( 1 2 ; - 21 i/2 - < 13 ) 

[47r 2 (n 2 + v 2 ) - a 2 ] ' 



Once again this expression requires the Fourier components b®(u,v) 
to be identically zero for u and v such that (u 2 + u 2 ) 1 / 2 < |a|/(2-7r). 
The large scale components of the field have infinite energy, and in 
this sense the general linear force-free boundary value problem has 
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no physical solution (Alissandrakis, 1981). In practice we use the Fast 
Fourier Transform (FFT) of the boundary field defined on an N by 
N grid over the region <x< 1, < y < 1. The smallest non-zero 
wavenumber that is determined from these 'observations' of the field 
(by analogy with magnetograph observations) is ii m i n = v m \ n = 
It follows that the energies determined from the discrete counterpart 
of (fTT?|) are finite for all N provided \a\ < a m ax = 2tt. 

Finally we note that the variation of (|13|) with a for small values of 
the force-free parameter is made obvious by a binomial expansion: 

-r== , ^ -r== , a 2 f f \b°Ju,v)\ 2 , , 

Ef(s)=Ef(0 ) + — 3 jj ^y^j^v 

3a 4 I I |!>!(u,i')l 2 , , 

+ ^JJ w^w iudv + (14) 

Hence for small a the non-potential component of the field scales with 
a 2 . 

2.3. Circuit model 

The current in a loop in the circuit model is taken to be the axial 
current contained within a radius R of the centre of a footpoint: 

[R- 

Ji = 2tt / J z r dr 
Jo 

« 2-nao 2 (\ - e -^/(2a?)^ ; (15) 

using the boundary conditions © or Q, and assuming that the foot- 
points are well separated. Similarly 

T 2 « 27raf3a 2 - e -« 2 /(2^I)^ ; (16) 

where (3 = B2/B1, and J z = [XqLJ z j B\. 

The energy of the single loop model in the half space ~z > is then 
estimated using (j3J): 

Ec = ^Ljl, (17) 

where 

Li = axCx, with C\ = In(8Si/H) - 7/4. (18) 
Similarly the energy of the two loop system can be expressed as 

E c = \Lj\ + ^L 2 7 2 2 + M12T1T2, (19) 
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where M± 2 is the mutual inductance. For the general case (loop major 
radii oi, 02, minor radii rq, f 2 , centres separated by di2, and axes ori- 
ented at an angle #12) the non-dimensional version of the interpolating 
formula Equation Q from Melrose (1997) may be used: 



_ m 2 1 al{c l c 2 ) l / 2 co S e 12 

M\2 = 



{a\ + ai) 2 + d 



12 



3/2 



(20) 



For the more restricted case of line-current loops with their axes aligned 
(812 = 0) the non-dimensional version of the exact Equation © may 
be used: 

Mi 2 = /i (aia 2 ) 1/2 

where 



2 

ki2 



ku K(k 



12 J 



-E(k 12 ) 



•■12 



(21) 



4aia 2 



"12 



2 

(ai + a 2 ) 2 + d 



(22) 



12 



3. Results 

3.1. One loop 

Figure 1 illustrates a single loop configuration, for the boundary condi- 
tions xi+ = (0.5, 0.1), xi_ = (0.5, -0.1), ai = 0.05, and 5=1. (These 
values are adopted as nominal values.) To construct this figure, the 
linear force-free field corresponding to (jllj) was calculated [using the 
FFT of Equation 0] on a 128 x 128 x 64 grid corresponding to the 
region 0<x<l,0<y<l,0<2<l. Field lines were traced for 
points at a radius cfi/2 = 0.025 from the positive footpoint. 

The free energy Epilx) — Ep(0) of the field can be estimated by 
taking the FFT of the boundary field with N gridpoints and then 
evaluating the integral (fT3|) using the trapezoidal rule. We label this 
estimate Wn- Figure 2 shows how the result depends on N, for the 
nominal choices of Xi±, (fi, and a. The upper panel plots Wjv versus N, 
for N = 32, 64, 128, 4096. The lower panel plots W N /W N/2 -l versus 
iV as a log-log plot. A condition for the sequence of values of Wn to con- 
verge is that Wn/Wn/2 approaches unity as iV — > 00. The lower panel 
suggests that this is occurring, and moreover Wn/Wn/2 — 1 ~ iV -1 . 
For the calculations in the paper we have used N = 512. The estimates 
of Wn for iV = 512 and iV = 4096 differ by about 0.5%, which suggests 
that iV = 512 is sufficiently accurate for present purposes. 
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The first question to be addressed is the variation of the energy 
estimates with I\ . To vary I\ we have varied a, while keeping a\ fixed 
[see Equation (|15|)]. We have also chosen R = a\, and used the nominal 
values of xi± and W\. The field estimate of the energy W512 and the cir- 
cuit estimate Eq were calculated for 50 values of a between 0.001a max 
and 0.99a max , with the values being equally spaced logarithmically. 
The results are plotted in Figure 3. The diamonds show the values of 
W512 and the crosses show the values of Eq. The vertical dashed line 
is a max and the horizontal dashed line is the estimate of Ep(0), the 
energy of the potential component of the field. The energy estimates 
for the non-potential component of the field and the circuit energy 
estimates are generally much smaller than the energy of the potential 
component of the field. The energy estimates for the non-potential 
component of the field exhibit a I x dependence (a 2 dependence) for 
most of the chosen range, but increase more rapidly for large I\ (a). 
This functional dependence is consistent with Equation (|14|) , The field 
and circuit energy estimates are within a factor of two of one another for 
most of the chosen range, but are out by more than a factor of 10 close 
to the maximum value of a. There is a simple physical interpretation 
for the observed dependence of the field energy on /. For small values 
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Figure 2. Upper: Variation of the free energy estimate Wn with N, for the single 
loop configuration. Lower: Convergence of the estimates. 



of the current the field associated with the current is small, and the 
magnetic field is dominated by the potential component. The current 
essentially follows the fieldlines of the potential field, i.e. there is a 
fixed geometrical path for the current. In this case the self- inductance 
defined by Equation (J2J) is independent of the current. For larger values 
of the current this is no longer the case: the field due to the current 
becomes comparable to the potential component of the field, and the 
path that the current takes then depends on the current. In particular, 
larger currents will tend to produce longer current paths, so the effective 
self-inductance will increase with current. In other words Equation (p|) 
increases with current. 

Another question is the dependence of energy on the geometry of 
the loop. The major radius a was varied between 0.05 and 0.25 for the 
nominal choices of the other parameters, and Figure 4 shows the results. 
The energy estimates W512 were divided by ^1 [ see Equation @] to 
produce effective self-inductances, which are shown by diamonds. Re- 
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Figure 3. Current dependence of energy estimates for the single loop configuration. 



suits are shown for a = 1.0 and a = 2.0 — the effective self-inductances 
for the larger value of a are slightly larger, for all values of major radius. 
The self-inductances L\ are indicated by crosses. This figure shows 
that the non-potential component of the field energy varies with major 
radius in a manner similar to that expected from the circuit model 
for the chosen range of values. The inferred self-inductance increases 
monotonically, as expected on the basis of the increase in area under 
the loop. However the precise functional dependence of the inferred self- 
inductance is somewhat different from that expected from the simple 
model L\. The similar results for the different values of a show that 
the inferred self-inductance is almost independent of the value of the 
current, as expected on the basis of the argument given above. 

The minor radius dependence of energy was also examined. Figure 5 
shows the effective self-inductance associated with the non-potential 
component of the magnetic field (diamonds) for a\ between 0.01 and 
0.1, and for the nominal values of the other parameters. The values 
of L\ are also plotted (crosses). This figure shows that the variation 
of the magnetic field energy with a\ is qualitatively consistent with 
the circuit model. The effective self- inductance decreases with a\, as 
expected from the decrease in area within within the loop as the minor 
radius of the loop increases. 
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Figure 4- Major radius dependence of self-inductance for the single loop configura- 
tion. 

3.2. Two LOOPS 

Figure 6 illustrates a particular two-loop configuration, for the bound- 
ary conditions x 1+ = (0.25, 0.575), xi_ = (0.25, 0.425), x 2+ = (0.825,0.5), 
X2_ = (0.675,0.5), ~&\ = (72 = 0.025, (3 = 1, and a = 1. These choices 
correspond to loops with major radii a\ = H2 = 0.075, with a distance 
d = 0.5 between their centres, and with axes at right angles to one 
another (9\2 = tt/2). 

There are many energy dependences which could be investigated for 
two- loop configurations. However, from the point of view of the circuit 
model the interest is with the energy associated with the interaction 
between the loops, described by the mutual inductance. The interaction 
energy is determined in the following way. The field energy VF512 for 
the two loops is calculated. Then the field energy ^512,1 for loop one 
is calculated (the free energy of the small scale linear force-free field 
with the boundary conditions of loop one). Similarly the field energy 
^512,2 for loop two is calculated. The energy W512 — W5i2,i — ^512,2 
is the interaction energy. Dividing this quantity by I1I2 defines an 
effective mutual inductance M\ 2 > which may be compared with circuit 
estimates. 

First we consider a comparison with the Stratton (1941) formula (|21j) . 
since it is an exact expression for mutual inductance (albeit for line 
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Figure 5. Minor radius dependence of self-inductance for the single loop configura- 
tion. 



currents). Recall that it is applicable only to the case of loops with 
aligned axes (612 = 0). The dependence of the mutual inductance 
on the major radius of loop two was investigated for loop parameters 
ai = 0.075, d = 0.5, a x = a 2 = 0.025, 12 = 0, /3 = 1, a = 1, with 
the major radius of loop two taking 20 values in the range 0.02 to 
0.2. Figure 7 shows the results. The points marked by diamonds are 
the mutual inductances calculated for the field following the procedure 
outlined above. The points marked with crosses are the mutual induc- 
tances given by the formula <|21JI . The inferred mutual inductances are 
qualitatively consistent with the circuit values, and a similar functional 
dependence is observed. 

Next we examine the Melrose (1997) interpolation formula 1)20(1 . 
which is more general than the Stratton (1941) formula, but is approx- 
imate. We considered a two-loop configuration with fixed parameters 
but with different values of #12, corresponding to a rotation of loop 
two with respect to loop one. Specifically the values 01=02 = 0.075, 
d = 0.5, ~&\ = ~&2 = 0.025, (3 = 1, and a = 1 were chosen, together 
with 21 values of 9i 2 in the range < #12 < vr. Figure 8 shows the 
results. The crosses are the circuit values. Following Equation ()20|) 
the circuit model exhibits a simple cosine dependence on 9i 2 , and in 
particular is zero for #12 = vr/2, because perpendicular loops have no 
flux linkage. The diamonds are the values of M^f determined following 
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Figure 6. A two-loop configuration with #12 = 7r/2. 

the procedure outlined above. The observed functional dependence is 
essentially cosine, and indeed the effective mutual inductance is very 
close to zero for 612 = 7r/2. However, the values of M 12 are considerably 
larger than the values given by the Melrose formula: the ratio of the 
values is close to eight over the range of 6*12 . This figure also shows the 
mutual inductance obtained using the Stratton (1941) formula for the 
special case 612 = as an asterisk. The Stratton value is slightly larger 
than M 12 for 9\2 = 0, and is considerably larger than the Melrose 
value. It is notable that the Melrose value differs from the Stratton 
value by more than an order of magnitude in this case. 



4. Discussion and conclusions 

In this paper calculated energies of non-potential components of linear 
force-free fields are compared with simple circuit model estimates, for 
boundary conditions suggestive of one and two magnetic loops. The mo- 
tivation is to test the circuit models, which are often invoked to estimate 
coronal magnetic field energies. In particular circuit models have been 
used to describe flare-related energy changes. In the Melrose (1997) 



ms.tex; 2/02/2008; 7:27; p. 12 



CIRCUIT ESTIMATES OF MAGNETIC FIELD ENERGIES 



13 



T3 
] 



0.08 



0.06 



0.04 



0.02 



0.00 




0.00 0.05 0.10 0.15 0.20 0.25 
Major radius of loop 2 (non — dim.) 



Figure 7. Variation of effective inductance with major radius of loop 2 for the 
two-loop configuration, and comparison with the Stratton (1941) formula. 



model, preferred magnetic configurations for flaring are identified on 
the basis of geometrical dependences of circuit energy estimates. 

In §3.1 results for a single loop configuration are presented. An 
important result is the a-dependence of the energy of the non-potential 
component of the field. For a range of small values of the force-free 
parameter a the free energy is found to be proportional to a 2 , or 
equivalently to be proportional to the square of the current in the loop, 
consistent with the circuit model. This result is explained mathemat- 
ically in terms of the a 2 dependence of the free energy appearing in 
Equation (|14|) for small a. For larger values of current (or a) the free 
energy increases more rapidly than the circuit prediction. This result 
is also expected for the two-loop configuration, since Equation (|14|) 
applies for any boundary conditions. A simple physical interpretation 
of this result is as follows. For small values of a the non-potential 
component of the field is small compared to the potential component, 
and the current essentially follows the fieldlines of the potential field. 
Hence as a is varied the geometry of the current does not change and 
the inductances [which are geometrical quantities — see Equations © 
and Q] are constant. This situation corresponds to the free energy 
scaling with the square of the current (or a 2 ). For larger values of a 
the non-potential component of the field becomes comparable to the 
potential component, and influences the path of the current. In this case 
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Figure 8. Variation of effective inductance with angle between loops for the two-loop 
configuration, and comparison with the Melrose (1997) formula. 



the geometry [and hence the inductances as defined by Equations © 
and depend on a, which corresponds to the free energy scaling with 
a higher power of a. Although demonstrated here for a class of linear 
force-free fields, this effect is expected to be quite general and to limit 
the accuracy of simple circuit models in application to any force-free 
field for large values of current. 

In §3.1 the dependence of the free energy of the field on the ge- 
ometry of the loop is also investigated (for a modest value of 75) . The 
functional dependences on minor radius and major radius are found to 
be qualitatively consistent with the circuit model. 

Section § 3.2 presents the results for the two- loop configuration. The 
interaction energy of the two loops is obtained by calculating the free 
energy of the loops together and then subtracting off the individual free 
energies of each loop. This procedure is used to calculate an effective 
mutual inductance M l2 for the field. This was then compared with sim- 
ple circuit model formulae, namely an exact expression for the mutual 
inductance of two axially-aligned line-current loops due to Stratton 
(1941), and an approximate expression for the mutual inductance of 
oblique loops with finite cross section due to Melrose (1997). The Strat- 
ton formula is found to qualitatively reproduce the observed variation 
of M 12 with the major radius of one loop. The Melrose formula is 
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compared with M 12 as one loop is rotated with respect to the other. 
The effective mutual inductance of the field reproduces the cosine de- 
pendence in the Melrose formula, and in particular the inferred mutual 
inductance is close to zero for the case of perpendicular loops. However, 
the calculated values of M^ 2 are almost an order of magnitude larger 
than the Melrose formula values. Comparison with the Stratton formula 
for the case of parallel loops confirms that, for the chosen parameters, 
the Melrose formula underestimates the actual interaction energy. The 
Melrose expression is an interpolation between known results, and so 
the discrepancy can be attributed to the approximate nature of the 
formula. 

In this paper linear force-free fields have been used because they 
are straightforward to calculate. It is likely that real coronal magnetic 
fields involve currents that are spatially highly concentrated, and hence 
require a nonlinear force-free model (or a description involving non-zero 
pressure forces) . In this case circuit models may provide quite accurate 
energy estimates because the configuration closely resembles a set of 
isolated circuits. However, it is also possible that in the nonlinear case 
the current path sensitively depends on the value of the current, in 
particular for large currents, and this may limit the accuracy of simple 
circuit models, as discussed above. 

Another point to note is that circuit models are really only useful 
if there are a small number of circuits which are easily identified. In 
this paper we have considered widely separated magnetic structures, 
which as a result have simple connectivity. In general coronal fields 
have a complicated topology, and it may be difficult to identify coronal 
current paths. 

Circuit models appear to provide reasonably accurate free energy 
estimates for coronal magnetic field configurations for a range of values 
of the current, and also qualitatively describe various functional de- 
pendences of the free energy (the dependence on loop major radii, the 
dependence on the angle between two loops, etc.). Simple circuit models 
are less accurate for very large values of the current, when the current 
path depends on the value of the current. Despite this reservation, the 
results of this paper confirm the general utility of circuit models, in 
particular in their application to the flare phenomenon. 
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